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LONG-RANGE PREDICTION OF FADING SIGNALS FOR WCDMA HIGH 
SPEED DOWNLINK PACKET ACCESS (HSDPA) 

CROSS REFERENCES TO RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional application Serial Number 
5 60/252,127 filed on November 20, 2000. 

FIELD OF THE INVENTION 

This invention relates to the field of wireless digital communications, and more 
particularly to digital signal processing for such signals. 

BACKGROUND OF THE INVENTION 

10 Wireless communications facilitates the delivery of information between the 

transmitter and the receiver without a physical wired connection. Such advantage 
translates to the freedom of mobility for the users and to the savings of wiring nuisance 
for the users. However, spectrum has become scarce resource as the usage of wireless 
communications for various applications becomes more popular. Therefore the 

15 efficiency of using spectrum presents challenges for the wireless industry. In order to 
maximize efficient spectrum utilization, various multiple access methods have been 
proposed to achieve the goal 

First generation cellular communications systems, Advanced Mobile Phone 
Services (AMPS) employed the Frequency Division Multiple Access (FDMA) method 

20 and provided voice communication services in the early days. Second generation cellular 
communications systems improved the spectrum efficiency by using more digital 
processing of signals and employed Time Division Multiple Access (TDMA) method in 
GSM and IS- 136 systems and Code Division Multiple Access (CDMA) method in IS-95 
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systems. While second generation systems typically provide two to five times voice 
capacity over the first generation systems, data capabilities of second-generation systems 
are very limited. 

A communication system where the transmitter has the side information feedback 
5 from receiver to transmitter was disclosed by Claude E. Shannon as early as in the 1950s. 
Channels with feedback from the receiving to the transmitting point are special case of a 
situation in which there is additional information available at the transmitter which may 
be used as an aid in the forward transmission system. Along with this line, a number of 
ideas have been presented which appeared to solve the problems in the fading channel. 
10 However, only recently the fading channel received a lot of attention due to the mobile 
wireless communications, particularly in the Code Division Multiple Access (CDMA) 
technology. 

SUMMARY OF THE INVENTION 

The present invention is an adaptive communication system, which supports 
15 higher peak data rate and throughput in digital wireless communications, compared with 
other non-adaptive systems. 

BRIEF DESCRIPTIONS OF THE DRAWINGS 

A more complete understanding of the present invention may be obtained from 
consideration of the following description in conjunction with the drawings in which: 
20 FIG. 1 is a high-level block diagram illustrating the principle of Long-Range 

Prediction and its application; 
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FIG. 2 is a diagrammatic representation showing the Channel State Information 
(CSI) obtained using either time-multiplexed pilot symbols (transmitted in DPCCH) or 
code-multiplex pilot channel signals (transmitted in CPICH); and, 

FIG. 3 shows a high-level block diagram of a WCDMA HSDPA system using 
5 Long-Range Prediction of Fast Flat Fading. 

DETAILED DESCRIPTION OF VARIOUS ILLUSTRATIVE EMBODIMENTS 

This invention digital is related to signal processing and system design, and more 
particularly to a mobile communication system for adaptive transmission in the radio 
frequency fading channel to improve the system capacity. The present invention is an 
10 adaptive system, which supports higher peak data rate and throughput in digital wireless 
communications, compared with other non-adaptive systems. 

One exemplary embodiment of the invention comprises three elements: the Long 
Term Prediction system for fast fading DS/CDMA mobile radio channel; the fast 
feedback system to enable the adaptive transmission; and new system blocks that are 
15 supported/enabled and changes in the existing 3GPP WCDMA system specifications. 

Fading of wireless signals is a deterministic process. One of the fundamental 
difficulties for the IS-95-B and IS-2000 standards lies in the fact that it is difficult for 
long duration of the frame structure to support fast channel information feedback. 

Principle of Long-Range Prediction 
20 In WCDMA, several adaptive transmission techniques, including adaptive 

modulation and coding, power/rate control, antenna diversity, ARQ, and others, are used 
for adaptation to rapidly time variant fading channel conditions. Since the channel 
changes rapidly, the transmitter and receiver are usually not designed optimally for 
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current channel conditions and thus fail to take advantage of the full potential of the 
wireless channel. By exploiting the time-varying nature of the wireless multi-path fading 
channel, all these adaptive schemes are trying to use power and spectrum more efficiently 
to realize higher bit-rate transmission without sacrificing the bit error rate (BER) 
5 performance. 

Referring to FIG. 1 there is shown a block diagram illustrating the principle of 
Long-Range Prediction and its application. Signal S(t) is coupled to a transmitter 102. 
The transmitter comprises an encoder 104, which is coupled to a modulator 106. The 
output of the transmitter 102 is X(t). Transmission channel 108 modifies the signal X(t) 

10 by multiplying the signal X(t) by the flat fading coefficient c(t) (as yet to be defined in 
Equation 1) and by the additive noise n(t), resulting in a modified signal 
y(t)=X(t) c(t)+n(t) which is detected by a receiver 110. The receiver 1 10 is comprised of 
a decoder 112 and a fading monitor & prediction using LRP section 114 which are 
coupled to the received signal y(t). The output of the fading monitor & prediction using 

15 LRP section 1 14 is coupled to the decoder 1 12 and a fast feedback channel 116 which is 
coupled to a modulation and coding selection (MCS) section 118. The output of the 
MCS section 1 18 is coupled to the encoder 104 and the modulator 106. 

Referring to FIG. 2 there is shown a diagrammatic representation showing the 
Channel State Information (CSI) 202 obtained using either time-multiplexed pilot 

20 symbols (transmitted in DPCCH) or code-multiplex pilot channel signals (transmitted in 
CPICH) 204. To implement the adaptive transmission methods, the channel state 
information (CSI) must be available at the transmitter. CSI can be estimated at the 
receiver and sent to the transmitter via a feedback channel. Feedback delay, overhead, 
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processing delay and etc are considered. For very slowly fading channels (pedestrian or 
low vehicle speed for most HSDPA applications), outdated CSI is sufficient for reliable 
adaptive system design. For faster speed, LRP is needed in order to realize the potential 
of adaptive transmission methods. These channel variations have to be reliably predicted 
5 at least several milliseconds (ms), or tens to hundreds of data symbols. Notice that one 
frame (15 slots) of WCDMA is 10 ms. The goal of LRP is to enable the adaptive 
transmission techniques. 

The present invention utilizes prediction of future fading conditions to improve 
the performance of WCDMA, especially for HSDPA applications. The present invention 

10 is a WCDMA system paradigm that uses the mechanisms of prediction of future fading 
conditions. The present invention is equally well suited for use with other system design 
such as CDMA2000. Of particular importance is how the new system paradigm 
improves the WCDMA system performance, especially high-speed packet access. 

Referring to FIG. 3 there is shown a high-level block diagram of a WCDMA 

15 HSDPA system using Long-Range Prediction (LRP) of Fast Flat Fading. In addition to 
the traditional system blocks, transmitter 302 and receiver 304 found in WCDMA 
HSPDA [3 GPP TR], new components including Fast Fading Monitor & Prediction Unit 
(FFMPU) 306, Reverse Link (RL) Fast Feedback Channel (RLFFC) 308, and Fading- 
Adaptive Unit (FAU) 3 10 are provided. 

20 The FFMPU 306 is simultaneously monitoring the current and predicting the 

future fast multipath fading using LRP. There are several LRP algorithms (to be 
discussed below) available for practical implementation. 
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The RLFFC 308 feedbacks some measured parameters describing the channel 
fading conditions from the mobile user equipment 304 (UE) to the base station 302 
(BTS). These parameters are measured in UE 304. 

The FAU 310 makes decisions on some selection on coding rate, modulation 
5 level, power allocation, multi-codes, number of rate matching bits required to fill a frame, 
ARQ, antenna diversity, scheduling, cell site selection, and etc. The FAU 310 can exist 
either in UE or BTS, depending on the final implementation complexity. 

The principle of FAU 3 10 is to adapt the selected system parameters to the rapidly 
changing fading channel conditions. The key feature of the system is the Long-Range 
10 Prediction ability of fading. Thus the transmitter 302 and receiver 304 have the accurate 
CSI parameters on future fading channel conditions by means of LRP. These CSI 
parameters include the maximum Doppler frequency shift. The availability of these 
forthcoming CSI parameters up to 15 slots/subframe in advance has made possible 
otherwise impossible new room in optimizing system design. Adaptation of the 
15 transmission parameters is based on the transmitter's perception of the channel conditions 
in the forthcoming time slots/subframes. Clearly, this estimation of future channel 
parameters can only be obtained by extrapolation of previous channel estimation called 
prediction. The channel characteristics have to be varying sufficiently slowly compared 
to the estimation interval. 
20 In the present invention, the inclusion of the LRP mechanism improves the 

WCDMA HSDPA system performance including supporting higher data rate. 

Adaptive Transmission Techniques used in Fading- Adaptive Unit (FAU) 
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The basic idea of adaptive modulation is to choose higher constellation size M of 
QAM (and therefore bit rate) for higher channel strength. Constant power and modulation 
size techniques suffer most BER degradation during deep fades. However fading channel 
spends most of the time outside deep fades. Thus adaptive modulation uses relatively 
5 high average constellation size (and bit rate) most of the time and avoid severe BER 
penalty by reducing the bit rate and using power efficient low modulation sizes (or 
turning off transmission entirely) during deep fades. The transmission load is shifted 
away from the deep fades and increases when the channel gets stronger. On the average, 
must faster bit rates relative non-adaptive techniques can be achieved without sacrificing 

1 0 the BER performance . 

The basic idea of adaptive channel coding is to select a code with lower rate when 
the channel is going into fade, and a higher rate when the channel becomes stronger. 
Punctured Turbo codes are used since they have superior performance and availability of 
a wide range of code rates without changing the basic structure of the encoder and 

15 decoder (codec). 

For adaptive transmitter diversity, the channel power of each transmitter antenna 
is monitored at the receiver, and the antenna with strongest power is selected. The 
diversity gain depends on how to accurately estimate the downlink propagation path 
conditions. LRP can improve this estimation. 

20 A critical fact for adaptive ARQ is that the transmission efficiency under flat 

Rayleigh fading conditions with smaller maximum Doppler frequency f d is higher than 
that AWGN channel conditions because long error-free length is more probable under 
fald Rayleigh fading conditions with smaller fd than under AWGN channel conditions 
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due to burstness of the error sequence. This is one of reasons that justify the use of ARQ 
or Hybrid ARQ in HSDPA. This fact also implies that "knowing" f d in advance of one 
future frame or future 10-15 slots/sub-frames, say, by means of LRP, seems to help the 
transmission efficiency using for a system using ARQ under flat Rayleigh fading channel 
5 conditions. When fd increases, transmission efficiency decreases because error-free 
length becomes short with increasing f d . Transmission efficiency depends on bit energy 
E b /N 0 . 

Scheduling of resources benefits from the knowing the future fading CSI and tries 
to avoid the transmission when channel is not in good conditions. The technique of the 

1 0 present invention will help reduce the scheduling delay and improve the throughput. 

Although space diversity is a very effective technique for compensating for rapid 
fading, it is helpless to compensate for log-normal fading or path loss due to distance. 
This requires so-called site diversity to obtain independent diversity paths by using plural 
base stations. In the case of Fast Cell selection, the UE selects the best cell every frame 

15 from which it wants to receive data on the HS-DSCH. HS-DSCH data is then transmitted 
to the UE from this cell only. UE can better select the best frame once UE knows the 
future fading CSI. 

If the fading CSI is known then the use of multi-code can be adaptively adjusted. 
Multiple Input and Multiple Output (MIMO) antennas seem to be sensitive to the 
20 fading CSI. The improved performance of LRP used for the fading CSI will definitely 
help MIMO antenna processing. 

LRP Algorithms in the FFMPU 
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LRP algorithms are known to those skilled in the art. A discussion of various 
algorithms can be found in LRP of Fading Signals by Alexandra Duel-Hallen, IEEE 
Signal Processing Magazine, May 2000, which is incorporated herein by reference as if 
set out in full. 

5 Signal Model 

Consider a low-pass complex model of the received signal at the user equipment 
r(t)=c(t) s(t) + I(t) Equation 1 . 

where c(t) is the flat fading coefficient (multiplicative), s(t) is the transmitted signal, and 
the I(t) includes the impact of the total interference resulting from the sum of M users , 
10 i.e. 

I(t) = Y j I ! (t) Equation 2. 

For the HSDPA case, we are interested in the downlink where the user equipment makes 
the measurement. I(t) can be modeled additive white Gaussian noise (AWGN). Let the 
transmitted signal at the base station be 

1 5 s(t) - J] b k " kT ) Equation 3. 

k 

where bk is the data sequence modulated using M-PSK or M-QAM, g(t) the BTS smitter 
pulse shape, and T the symbol delay. At the output of the matched filter and sampler, the 
discrete-time system model is given by rk = bk Ck + z^ where Ck is the fading signal c(t) 
sampled at the symbol rate, and Zk is the discrete AWGN process I(t). In general, the 
20 sampling rate represented by subscript n differs from the data rate represented by k 
throughout this paper. Usually, c(t) and Ck can be modeled as a correlated complex 
Gaussian random processes with Rayleigh distributed amplitudes and uniform phases. 
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Using the pilot-aided signals in WCDMA, the receiver can correctly detect the symbol bk 
. Then by multiplying the received samples by the conjugate of bk , the modulation can be 
removed, yielding 

ik = Ck + zjk Equation 4. 

5 where zik is still an AWGN with the same variance as Zk. 

The derivation of this prediction method is based on a physical description of the 
fading signal. In this section, the mathematical description of the interference pattern 
from the point of view of the mobile is primarily considered. The fading coefficient at the 
receiver is given by a sum of N Doppler shifted signals 

10 c(0 = XV WM Equation 5. 

n=\ 

where (for the n-th scatterer) A n is the amplitude, f n is the Doppler frequency, and <j) n is 

the phase. The Doppler frequency is given by 

f n = fc (v/c) cos(ocn) Equation 6. 

where f c is the carrier frequency, v is the speed of mobile, c is the speed of light, and a n is 
1 5 the incident angle relative to the mobile's direction. Due to multiple scatterers, the fading 

signal varies rapidly for large vehicle speeds and undergoes "deep fades". 

The fading signal Ck in Equation 4 is predicted by decomposing it in terms of the 

N scattered components. If the parameters A n , f n > and a n in Equation 5 for each of the 

scatterers were known and remained constant, the signal could be predicted indefinitely. 
20 In practice, they vary slowly and are not known a priori. Assume that the propagation 

characteristics will not change significantly during any given data block. Therefore, these 

parameters are modeled as approximately constant or change slowly varying for the 
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duration of the data block. To predict the fading signals, spectral estimation followed by 
linear prediction and interpolation is employed. Estimation of the power spectral density 
of discretely sampled deterministic and stochastic processes is usually based on 
procedures employing the Discrete Fourier Transform (DFT). Although this technique 
5 for spectral estimation is computationally efficient, there are some performance 
limitations of this approach. The most important limitation is that of frequency 
resolution. The frequency resolution Af=l/f s of the N-point DFT algorithm, where f s is 
the sampling frequency, limits the accuracy of estimated parameters. These performance 
limitations cause problems especially when analyzing short data records. 
10 Many alternative Spectral Estimation Techniques have been proposed within the 

last three decades in an attempt to alleviate the inherent limitations of the DFT technique. 
What follows are several practical alternative embodiments, considering the specific 
application to HSDPA. 

Maximum Entropy Method (MEM) 

15 The Maximum Entropy Method (MEM) for the prediction of the fast fading 

signal, is also known as the All-poles Model or the Autoregressive (AR) Model and is 
widely used for spectral estimation. The reason this technique was chosen is that the 
MEM has very nice advantage of fitting sharp spectral features as in the fading channel 
due to scatterers (Equation 5). Furthermore, MEM is closely tied to Linear Prediction 

20 (LP), which is used to predict future channel coefficients. Using MEM, the frequency 
response of the channel is modeled as: 

H{z) ~ Equation 7. 

7=1 
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This model is obtained based on a block of samples of the fading process. Note that the 
samples have to be taken at least at the Nyquist rate, which is twice the maximum 
Doppler frequency, fd. Moreover, the accuracy of the model depends on the number of 
samples in the given block. The dj coefficients are calculated from the poles of the power 
5 spectral density. The dj coefficients in Equation 7 are also the linear prediction 
coefficients. The estimates of the future samples of the fading channel can be determined 
as: 

p 

c n = ^]d J c n _j Equation 8 , 

Thus, c n is a linear combination of the values of c n over the interval [n-p, n-1]. Since 
10 actual channel coefficients are not available beyond the observation interval, earlier 
sampling estimates c n ^ } , can be used instead of the actual values c n .j in Equation 8 to 

form future estimates c n , or the samples can be updated adaptively below. 

Note that the channel sampling rate utilized for LP is much lower than the symbol 
rate, 1/T. Therefore, to predict the fading coefficients, Ck in Equation 4, associated with 
15 transmitted symbols, interpolation is employed as discussed. In this interpolation 
process, four consecutively predicted channel coefficients are interpolated by a Raised 
Cosine (RC) filter to generate estimates of the fading coefficients, c k , between two 

adjacent predicted samples at the data rate. 

Interpolation is preferred to oversampling of the fading channel to obtain the 
20 fading coefficients at the data rate. If oversampling is employed, MEM will require a 
larger number of poles and consequently the complexity will increase. 
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The prediction method can be combined with tracking and transmitter signal 
power adjustment. The channel samples taken during the observation interval are sent to 
the transmitter, which applies linear prediction to compute the coefficients and 
interpolates to produce predicted fading values at the data rate. Note that this feedback is 
5 not going to introduce significant delay since the sampling rate is much lower than the 
data rate. Then, the transmitter sends the data bits, bk, by multiplying them with the 
inverse of the c k values. While this is not the optimal method for transmission over the 
time varying channel, it still achieves significant gains relative to the case when power 
compensation is not employed at the transmitter. At the output of the matched filter and 
10 sampler, the new modified discrete-time received signal is given by 

yk=T Lb k +z k Equation 9. 

Ck 

where z k is discrete-time AWGN. Define a k = — . As the prediction gets better, the 

c k 

value of ak goes to 1. When ak=l, i.e., perfect estimation, our fast fading channel becomes 
the AWGN channel. 

15 The Least Mean Squares (LMS) adaptive algorithm is employed to track the 

variations in a k . Given the received signal (Equation 9), the LMS algorithm is performed 
at the data rate as 

^jt+i = %k + A (y k -VkY Equation 1 0. 

where \x is the step size, y k =a k bk . This tracking is employed to perform coherent 
20 detection at the receiver, as well as to update the estimate of the fading at the sampling 
rate. The new fading sample is computed as c k = a k c h and send back to the transmitter at 

-13- 


4648-105 US 

the sampling rate. The transmitter uses this updated estimates in (8) to predict future 
fading values, rather than relying on previous estimates. This adaptive algorithm enables 
us to reduce the prediction error described earlier and to approximate the performance of 
the AWGN channel. 
5 Root-MUSIC 

Root-MUSIC is especially useful, in that it has two desirable features: high 
resolution and no need for spectral peak finding. 

A K-by-K sample correlation matrix can be constructed from output data in 
Equation 4, i.e., 

10 R=G*G H Equation 11. 

Where G is the forward-backward data matrix constructed from output data in Equation 
4. Assuming the number of sinusoids P (typically P=8-10) is known, then the noise 
subspace is obtained as Span{V n }, =[ V p +i V p+ 2 V K ] Where V n consists of the K-P 
smallest eigenvectors of R. Let Q= V n V n H and 

K-i K-i 

15 c,=X<&,* + ' and <L,=£& + , fjfe M / = 0,1,2,...,* -1 Note that cVc, ; 

and forms the polynomial equation c_ K+] + c_ K+2 z' ] + ... + c 0 z~ K+] + ... + c K ^z~ 2{K ~ x) - 0 . 
Solving this equation gives 2(K-1) roots having reciprocal symmetry with respect to the 
unit circle. Denote the P roots that outside and also nearest to the unit circle as zi, Z2, 
....z p . Then the frequency estimates are given by fi=arg(zi)/27i ? i=l,2,...,P, where arg(z,) 
20 denotes the principal argument (in radians) of z,. Root-MUSIC needs to know the number 
of the sinusoids a priori. So-called root location constraints can be used to avoid this 
problem. 
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Once the frequency estimates have been obtained, the complex amplitudes 
E^Aie* can be found by linear least square (LS) fit of the following matrix-vector 
equation A E= [ai a 2 ...a p ] E=g, where ai=[l e* 2%n ... e? 2nm ] T for i=l,2,...,P, E=[Ei 
E 2 ...E P ] T is the complex vector to be found, and g=[g(0) g(l) ...g(N)] T . The LS solution 

5 of the above equation is given by E =A # g, where A # =(A H A) _i A H is the pseudo inverse of 
A. In this way, the parametric sinusoidal model for the fading process is obtained. 
Fading prediction can be done by this method. 

MMSE AR Method 
MMSE prediction of the flat fading channel is used for the AR model 
10 ESPRIT 

Performance Bounds 
The performance of the method is described as following. For ak=l , i.e., perfect 
estimation, the average probability of bit error for BPSK is given by 

P e = Q(^2n) Equation 12. 


15 where yb is the SNR and Q(x) is defined as Q(x) = -j= \e ' 2/2 at . Since this performance 


is achieved with perfect prediction and this curve forms the lower bound for our system. 
If there is no correction at the transmitter, the received signal is given by Eq. (4). Since Ck 
is approximately Rayleigh, the average probability of bit error for the Rayleigh fading 
channel is found as 


20 p=i(i Equation 13. 
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Equation 13 forms the upper bound of the proposed method. The expected realistic 
performance should lie between the upper bound and lower bound. 

For the QAM similar curves are obtained. For square-QAM, carrier regeneration 
using pilot-aided signal is essential. Gray encoding with absolute phase coherent 
detection can be applied. The BER for Gray-encoded 16QAM and 64QAM is, 
respectively, for AWGN given by 


Equation 14. 


For Rayleigh fading channel, it is seen that 


P = 1 

1 e\60AM 8 


1-- 


1 + 


2n 


p =^ 

r e640AM 24 


1-- 


1 + - 


Equation 15. 


10 Numerous modifications and alternative embodiments of the invention will be 

apparent to those skilled in the art in view of the foregoing description. Accordingly, this 
description is to be construed as illustrative only and is for the purpose of teaching those 
skilled in the art the best mode of carrying out the invention. Details of the structure may 
be varied substantially without departing from the spirit of the invention and the 

1 5 exclusive use of all modifications, which come within the scope of the appended claim, is 
reserved. 
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For example, although the inventive concept was illustrated herein as being 
implemented with discrete functional building blocks, the functions of any one or more of 
those building blocks can be carried out using one or more appropriately programmed 
processors, e.g., a digital signal processor. It should be noted that the inventive concept 
5 is equally well suited for other wireless systems. 

In one exemplary embodiment, the present invention supports higher peak data 
rate and throughput, compared with other non-adaptive systems. In yet another 
exemplary embodiment, the present invention can be supported by the existing 3 GPP 
WCDMA system structure, particularly the frame/slot structure. The present invention is 
10 equally valid for use with other similar systems where the frame structure supports the 
fast feedback from receiver to transmitter point. Once the principle of fading adaptation 
is established, each related part of the mobile communications system can be improved. 

While various terms and abbreviations are defined in this application, and would 
be clearly understood to and understood by one skilled in the art, attention is drawn to the 
15 above referenced publications for further details and descriptions. 
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